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ABSTRACT 

We present results from a low-resolution spectroscopic survey for 21 galaxy clusters at 0A < z < 

0.8 selected from the ESO Distant Cluster Survey (EDisCS). We measured spectra using the Low- 
Dispersion Prism (LDP) in IMACS on the Magellan Baade telescope and calculate redshifts with a 
precision of a z = 0.006. We find 1,602 galaxies that are brighter than R = 22.6 in the large-scale 
cluster environs. We identify the galaxies expected to be accreted by the clusters as they evolve to 
z = 0 using spherical infall models, and find that ~ 30-70% of the z = 0 cluster population lies outside 
the virial radius at z ~ 0.6. For analogous clusters at z = 0, we calculate that the ratio of galaxies 
that have fallen into the clusters since z ~ 0.6 to that which were already in the core at that redshift 
is typically between ~ 0.3 and 1.5. This wide range of ratios is due to intrinsic scatter and is not 
a function of velocity dispersion, so a variety of infall histories is to be expected for clusters with 
current velocity dispersions of 300 <cr< 1200 km s _1 . Within the infall regions of z ~ 0.6 clusters, 
we find a larger red fraction of galaxies than in the field and greater clustering among red galaxies 
than blue. We interpret these findings as evidence of “preprocessing”, where galaxies in denser local 
environments have their star formation rates affected prior to their aggregation into massive clusters, 
although the possibility of backsplash galaxies complicate the interpretation. 


1. INTRODUCTION 


Although a relationship exists between the evolu¬ 
tion of galaxies and their environment, as demonstrated 
by correlat i ons be tween density and galaxy colo r (e.g. , 


Hogg et al.1 20041), star-formation (1Lewis et al. 


2002, 

Gomez et al.1 20031), and morphology (Dressier 

1980; 

Postman & Gelled 1984|), the physical processes 

that 


drive these changes and the connection between those 
processes and environment are not established. While 
the cores of clusters are the final resting place for quies¬ 
cent galaxies and are where these trends were discovered, 
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the key to understanding the implicit quenching of star 
formation and morphological transformation is to study 
galaxies in the environment where they are being trans¬ 
formed, not where they ultimately reside. 

Quenching and morphological transformation do not 
occur primarily in the cores of clusters, at least not at 
redshifts < 1. The de crease in star formation sets in 
at se veral virial radii (|Lewis et all 120021 I Gomez et al.1 
120031 b and the increase i n the SO fraction s ince z ~ 
0.5 ( Dressier et al.1 119971 : iFasano e t all 120 00(1 is most 
dramatic in less -massive clusters (Poggi anti et all [20091 
I Just et all I2010T) . Environmentally-driven evolution oc¬ 
curs primarily at intermediate densities, which should in¬ 
clude the environs outside the cluster virial radius. Such 
effects are predicted in si mulations out to as many as 
five virial radii dBahe et al.ll2013l l. Establishing the size 
and characteristics of the infalling galaxy population will 
therefore constrain the path to transformation. 

As a result of this line of thought, a number of stud¬ 
ies have begun to target the outskirts of z> 0.5 mas¬ 
sive clusters fe.g.. iMoran et al.1 [20071 iPatel et all 120111 : 
IQemler et al.1120131b However, such studies have been 
limited to a few clusters (~ 10), making general conclu¬ 
sions difficult to reach given the variation in properties 
from cluster to cluster. Because of the high masses of 
these targeted clusters, and correspondingly large virial 
radii, some of these studies do not probe very far past the 
virial radius and may miss a significant fraction of the in¬ 
falling galaxies. Furthermore, such clusters are also rare, 
hence the infalling population of more typical clusters 
has not been explored. This bias may lead to an incom¬ 
plete pict ure, given the cluster - mass dependenc e of SO 
evolution ([Poggianti et al.ll2009HJust et al.1 1 20101) . 

A fundamental difficulty in studying cluster infalling 
populations is the contamination of interloping fore¬ 
ground/background galaxies, an issue that becomes more 
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important at larger cluster-centric radii where the rela¬ 
tive fraction of interlopers is larger. The studies listed 
above use spectroscopic redshifts for this purpose, but 
this approach requires significant telescope time and is 
thus limited to those few clusters. The alternative ap¬ 
proac h using photometric redshifts (e.g.. [Kodama et ahl 
12001 comes with much lower observational cost, but 
photometric redshifts are insufficiently precise to se¬ 
curely associate a galaxy to a particular cluster, where 
cdz ss 500 km s -1 resolution is needed. 

We adopt a hybrid approach. We isolate the infalling 
galaxy population of 21 clusters at 0.4 < z < 0.8 
using the Low-Dispersion Prism (LDF0) installed in 
the Ina mori-Magellan Areal Camera and Spectr ograph 
('IMACS JBigelow et al.lll998tlDressler et al.ll2006l f on the 
6.5m Magellan Baade telescope. With these data, we 
measure the number of galaxies these clusters will ac¬ 
crete by z = 0 to establish how many galaxies may be 
influenced by the accretion process. We also measure the 
scatter in this number to estimate the range i n accre¬ 
tion h istories. We compare models ('e.g.. iPoggianti et al.l 
I2006T > that predict the amount of mass accreted by these 
clusters to our observations. Finally, we measure the 
optical properties and clustering amplitude of infalling 
galaxies to quantify the amount of evolution that takes 
place outside the virial radius. 

This paper is organized as follows. In §2 we describe 
our sample selection and the sample’s basic properties, 
and in §3 and §4 we present the imaging and spectro¬ 
scopic data, respectively. In §5 we analyze our clus¬ 
ters using mass infall models and quantify the num¬ 
ber of galaxies and optical properties of the infalling 
population. We conclude in §6. All magnitudes in 
this paper are in the AB-system; to convert these to 
the Vega system, subtract 0.02, 0.06, 0.23, 0.45, and 
0.55 from the AB magnitudes for the BVRIz bands, 
respectively. Throughout the paper, we adopt Hq = 
70 km s -1 Mpc -1 , Do = 0.3, and Da = 0.7, and all 
cosmology-dependent quantities taken from other studies 
also use these values. We approximate the virial radii of 
our clusters as J? 200 > the radius inside which the enclosed 
density is 200 times the critical density of the Universe 
at that redshift. 


2. SAMPLE 

Our sample consists of 21 galaxy clusters. We include 
16 of the 20 gala xy c l usters in the ESO Distant Cluster 
Survey fEDisCS: IWhite et alj|2005h : see 44.11 for details 
of the four clusters not observed. We also include the 
seven clusters found serendipitously in this survey. Of 
these 23 LDP-observed clusters, two are removed from 
the analysis for reasons given in 44.21 resulting in 21 clus¬ 
ters in our final sample. We present basic information on 
the clusters in Table 1. 

The EDisCS clusters were drawn from candi- 
dates in the Las C ampanas Distant Cluster Survey 
(jGonzalez et al.ll2001h identified as surface brightness en¬ 
hancements in the image background. They lie in a band 
from ss 10-14 hr in right ascension and « —13 to —11 de¬ 
grees in declination. They span a redshift range from 
2 = 0.4 to 0.8 and cover a spread in velocity disper- 

1 5 Designed by S. Buries for use by the PRIMUS redshift survey 
ICoil et al .1120 1 if ). 


sion ( a) in the range sa 200-1200 k m s -1 (jHallidav et al.1 
12004 IMilvang- Jensen et al.l 120081) . a wider range of a 
than other cluster samples at these redshifts and more 
representative of the prog enitors of 2 ~ 0 clusters 
(jMilvang-Jensen et al.|[2008l) . 

We have a variety of data on the cluster cores (the 
central ~ 6.5' x 6.5' FOV), with deep (/< 25) opti¬ 
cal imaging from the FOcal Reducer and low disper¬ 
sion Sp ectrograph (FORS2) on the Very Large Telescope 
(VLT) (1 White et al.1 120051) . near-infrared imaging from 
the So n OF ISAAC (SOF T) at the New Technology Tele¬ 
scope dWhite et al.ll2005R . and optical VLT spect roscopy 
(|Hallidav et al.ll2004 IMilvang- Jensen et al .1120081) . weak- 
lensi ng maps (fClowe et aT 1 20061), gala xy morpholo¬ 
gies (|Desai et al.ll2007t~ Simard et al Il2009h . fundamental 
plane parameters (ISagli a et al.1 201011 . bri ghtest cluster 
galaxy identifications ( While ve t al.l 20 08h and MIPS - 
based star-formation rates ( Finn et al.l 120051 [2010). 
Wide-field imaging in the mid-infrared with MIPS (~ 
50' x 20' FOV) and ultra-violet with GALEX (ss 38' 
radius FOV) also exist for cluster subsets, but do not 
appear in this study. 


3. WIDE-FIELD IMAGING DATA 

We use wide-held (~ 30' x 30') imaging of our clus¬ 
ters to identify targets for our LDP masks and to mea¬ 
sure galaxy magnitudes and colors, which are used for 
the redshift-htting portion of the LDP pipeline as well 
as for characterizing the galaxies. Our photometry is 
measured from BVRIz images, with VRI data from 
the Wide Field Imager (WFI) instrument on the 2.2m 
Max Planck Gesellschaft/E uropean Southern Observa- 
tory (MPG/ESO) telescope (jBaade et al.lll999l ) and Bz 
data from MOSAIC on the Cerro Tololo Inter-American 
Observatory (CTIO) Blanco telescope, which have 34' x 
33' and 36' x 36' FOVs, respectively. Not all clusters have 
been observed in all 5 bands. Our entire sample has VRI 
data, while some clusters appearing in IGuennou et al.l 
(|2010ll have either B or 2 , or both (see Column 9 of Ta¬ 
ble 1). Details on the imaging data are given below. 


3.1. VR-, and I-band Data from WFI 

We redu ce the raw images u sing the tec hniques de¬ 
scribed by iClowe fc Schneideil (|200ll l2002lf . which in¬ 
volve bias subtracting and flat fielding each chip sepa¬ 
rately and removing fringing in the R- and /-band im¬ 
ages. We calculate astrometric solutions for the images 
by comparing the image centroids of U.S. Naval Obser¬ 
vatory (USNO) reference stars, and use the utility imwes 
to write a new world coordinate system (WCS) header 
based on those matchesFl This procedure results in an 
RMS position per star of ~ 0'/3 relative to the USNO co¬ 
ordinates. For C11354.2-1230, this method failed to con¬ 
verge , so we define the astrometry using SCAMP dBertinl 
l2006f h The astrometric precision is « 0"5 for this field. 

We create photometr ic catalogs using SExtractor 
(jBertin fe Arnoutsl 1199611 . We detect sources in the 
seeing-matched A-band image, requiring at least 12 ad¬ 
jacent pixels containing flux > 5 ctrms above the back¬ 
ground. Photometry is performed in two-image mode 

16 Originally written at the University of Iowa, but since adapted 
and amplified by Doug Mink at the Smithsonian Astrophysical 
Observatory (http://tdc-www.harvard.edu/wcstools/imwcs/). 
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TABLE 1 

LDP-Observed EDisCS Clusters 


Field 

Cluster ID 

RA 

Dec 

z 

<7 

R200 

M200 

Imaging 

Seeing (") 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

1 

C11018.8-1211 

10:18:46.8 

-12:11:53 

0.4734 

486+gf 

0.931?;)) 

1.5318:81(14) 

VRI 

1.20 

2 

C11037.9—1243 

10:37:51.2 

-12:43:27 

0.5783 

319tg 

0-58l°;)° 

4-06l?;||(13) 

BVRIz 

2.10 

3 

C11037.9—1243a 

10:37:52.3 

-12:44:49 

0.4252 

537!« 

1 QC + 0.09 
l.UD—0.09 

2.1218;?? (14) 

BVRlz 

2.10 

4 

C11040.7-1155 

10:40:40.4 

-11:56:04 

0.7043 

418 !$g 

0.70lg;gl 

8-47ll;|8 (13) 

BVRIz 

1.45 

5 

C11054.4—1146 

10:54:24.5 

-11:46:20 

0.6972 

589!)? 

0.991°;)? 

2.38l);°?(14) 

BVRIz 

1.20 

6 

C11054.7—1245 

10:54:43.6 

-12:45:52 

0.7498 

504l“ 3 

0.82!°;)? 

1.441);?) (14) 

BVRIz 

1.25 

7 

C11059.2-1253 

10:59:07.1 

-12:53:15 

0.4564 

510 !?? 

0.9918;)? 

1.7818;?? (14) 

VRI 

1.05 

8 

C11103.7—1245a 

11:03:34.9 

-12:46:46 

0.6261 

336tg 

0 15Q+ 006 

u.oy_ 0- 07 

4.611)21(13) 

BVRI 

1.15 

9 

C11103.7—1245b 

11:03:36.5 

-12:44:22 

0.7031 

252tg 

0 42 +0 11 
u -^ z —0.14 

1.86l);||(13) 

BVRI 

1.15 

10 

C11138.2—1133 

11:38:10.3 

-11:33:38 

0.4796 

732tS 

1 40+ 014 
1 - 4U _0.15 

5.20l);«° (14) 

BVRI 

1.15 

11 

C11216.8-1201 

12:16:45.1 

-12:01:18 

0.7943 

I— 1 

o 

00 

1 + 
-j -i 

-4 CO 

1 fi1+ 012 

1-O1-0.12 

i-i6l8;II(i5) 

BVRI 

1.20 

12 

C11227.9—1138 

12:27:58.9 

-11:35:13 

0.6357 

574 !)? 

1 on +013 

1 - uu -0.13 

2.2918;?) (14) 

BVRI 

1.25 

13 

C11227.9—1138a 

12:27:52.1 

-11:39:59 

0.5826 

341 ±21 

0 fii+ 0 08 
u,D1 -0.08 

4.95l?°°(13) 

BVRI 

1.25 

14 

C11232.5—1250 

12:32:30.5 

-12:50:36 

0.5414 

1080l“ 9 

1 qq+0.22 

i-yy_o.i6 

1.6ll°;®?(15) 

VRIz 

1.05 

15 

C11301.7-1139 

13:01:40.1 

-11:39:23 

0.4828 

6871“ 

1 *31 “I -9 - 

J--oi_o.i6 

4.29l);7?(14) 

VRI 

1.15 

16 

C11301.7—1139a 

13:01:35.1 

-11:38:36 

0.3969 

39ltfjj 

0.78l°;) 3 

8.32l);™(13) 

VRI 

1.15 

17 

CU353.0-1137 

13:53:01.7 

-11:37:28 

0.5882 

6661) 3 ® 

1-I9l?;?5 

3.67!?;)? (14) 

VRI 

1.20 

18 

C11354.2—1230 

13:54:09.7 

-12:31:01 

0.7620 

6481)?® 

1,uo -0.18 

3.05!);)? (14) 

BVRIz 

1.66 

19 

C11354.2—1230a 

13:54:11.4 

-12:30:45 

0.5952 

4 QQ+ 95 

4 * 5 * 5 — 104 

0 77+ 017 
u - ‘ * -0.19 

i.oo!8;?8(i4) 

BVRIz 

1.66 

20 

C11411.1-1148 

14:11:04.6 

-11:48:29 

0.5195 

7101)11 

1 qq + 0.23 
i.oo_Q 25 

4.63!?;?° (14) 

VRI 

1.45 

21 

C11420.3-1236 

14:20:20.0 

-12:36:30 

0.4962 

2181$ 

0 41 +0-08 

u, ^ :1 -0.09 

1.36lg;?)(13) 

VRI 

1.00 

22 

C11103.7-1245 

11:03:43.4 

-12:45:34 

0.9586 

534+ 101 
00 -120 

0 77+ 0 - 15 
u -' '-0.17 

1.52l)'°)(14) 

BVRI 

1.15 

23 

C11138.2—1133a 

11:38:08.6 

-11:36:55 

0.4548 

5421?? 

1 O'S - * -012 

l.UD-o 14 

2.14l°;??(14) 

BVRI 

1.15 


Note. — (1) cluster field; (2) cluster name; (3,4) J2000 right ascension (hours) and declination (deg); (5) cluster redshift 
(Halliday ct al. 2004; Milvang-Jensen et al. 2008); (6) cluster velocity dispersion in km s~ 1 (Halliday ct al. 2004; Milvang-Jensen et al. 
2008); (7) cluster virial radius in Mpc; (8) cluster virial mass in units of Mq with power of ten in parentheses (using Equation 10 
of lFinn et al.l[2005fl : (9) wide-field imaging bands observed in each field; (10) the effective WFI seeing after smoothing the images to 
match the band with the poorest seeing for that cluster 

for the other bands. Given the wide FOV, we correct for 
Galactic extinction differentially across the field. The 
color excess, which is directly proportional to the extinc¬ 
tion, across a given field varies by ft 0.01-0.02. We de¬ 
termine E(B — V) at each p h otom etric source using the 
dust maps of ISchleg el et al.| (|1998D. and interpolate the 
extinction curve of lCardelli et al.1 (|1989D to the effective 
wavelength of each bandpass to determine the extinction, 
assuming Ry = 3.1. 

To match the point-spread functions (PSFs) among 
bands so that aperture-matched magnitudes probe the 
same region of the galaxy, we smooth the images with 
a gaussian kernel selected to match the image with the 
largest seeing for that field (often the V or I band). The 
resultant effective seeing is typically 1"2 (full-width half 
max, or FWHM) for the different fields, except C11037.9- 
1243, which has seeing ft: 2". For most clusters, the im¬ 
age quality or effective seeing varies by less than O'/l (< 

0.5 pixel) over the image; for C11227.9-1138, C11232.5- 
1250, C11353.0-1137, CU354.2-1230, and CU411.1-1148 
it varies < 0'.'2 (< 1 pixel). 

The WFI data were taken under non-photometric con¬ 
ditions, and therefore are poorly calibrated. We adjusted 
the photometric zero points (ZPs) for the VRI data us¬ 
ing well-calibrated and deep VLT images taken as part of 
the original EDisCS program. To determine the ZPs, we 
first cross-correlate stars from the WFI images with those 


from the VLT using a 0"5 matching threshold, resulting 
in ~ 20-100 matches per field. We compare the non- 
extinction-corrected VLT magnitudes of these matches 
with their counts in S^-radius apertures on the WFI im¬ 
ages, and use linear regression to calculate color terms 
(a\) of the following form for each of the bands: 

ZPw fi,a = ■Z’^Vlt.a + a\(V — R) wfi,a- (1) 

A first guess for the V — R color yields WFI ZPs with 
which we calculate new VRI magnitudes using 

wwfi,a = -2.51og 10 (countswFi,A) + ZR WF i,a- (2) 

These in turn give new V — R colors. This process is 
iterated until convergence. The uncertainties in VRI 
ZPs are 0.12, 0.07, and 0.12, respectively. 

R-band imaging from the VLT was not available for 
C11018.8-1211. CU059.2-1253, C11232.5-1250, CU301.7- 
1139, C11353.0-1137, C11411.1-1148, and CU420.3-1236. 
In what directly follows, all bandpasses refer to the VLT 
filters. We estimate R-band magnitudes from synthetic 
R-band magnitudes that we obtain by fitting the BVIK 
spectral ene rgy distributions (SEDs) with stellar tem¬ 
plates from iHauschildt et al.1 (119991 ). We also use this 
methodology on clusters with R-band photometry to as¬ 
sess its accuracy, finding that the absolute value of the 
difference between predicted and observed mean R-band 
magnitude is < 0.02. For the seven clusters without VLT 
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R-band imaging, we use the R-band magnitudes derived 
in this way for the iterated scheme described above when 
measuring ZPs. 

To calculate colors, we use fixed apertures of l"-radii 
to maximize signal-to-noise; using a larger aperture in¬ 
troduces more noise into our color measurement. We 
estimate the total magnitude using the FLUX_AUTO 
measurement from SExtractor, which fi ts so urces with 
an ellipse following the method of iKro itfUIO). 

We estimate photometric errors by placing 10 3 back¬ 
ground apertures on each image with radii ranging from 
1 to 6 pixels (« 0.2-3"), avoiding sources by using the 
SEGMENTATION output of SExtractor. We fit the 
RMS fluctuation of counts in each aperture as a func¬ 
tion of aperture size to determine the error at 1" and, 
for the AUTO magnitudes, at the Kron radius. These 
errors are typically « 0.01-0.02 for most of the galax¬ 
ies that appear in this paper, although at the magnitude 
limit (see below), the V-, R-, and J-band errors approach 
0.03, 0.02, and 0.03, respectively. Errors on V — I color 
are « 0.1 or less. 

Comparison with the VLT magnitudes suggests an 
RMS precision that varies with limiting magnitude, in 
that the RMS increases at faint magnitudes. If we only 
include WFI galaxies with Rauto < 23.3 (our photomet¬ 
ric completeness limit; see below), the RMS precision is 
0.12, 0.07, and 0.12 for the V, R and I total magnitudes, 
respectively. This includes the Poisson uncertainty of 
counts in the aperture. 

3.2. B- and z- Data from MOSAIC 

We use B- and g-band d ata for nine of our clusters ob¬ 
tained bv iGuennou et all (12010H with the CTIO Blanco 
telescope using MOSAIC. These data were r educed w ith 
the MIDAS, SCAMP and SWarp pack ages (jBanse et al.1 
119881 IBertin et al.1120021 IBertinl 1200611 . Exposure times 
for the B and z data a r e 11 x 600s and 18 x 800s, respec¬ 
tively. IGuennou et al.1 (]2010H describe the data in more 
detail. 

The B- a nd z-band ZPs have e rrors of 0.09 and 0.07, re¬ 
spectively (IGuennou et al.ll2010ll . The B-band ZPs were 
corrected for galactic exti nctio n using a single E{B — V) 
value per field from the lSchlegel et al.1 (119981) dust maps; 
we “de-correct” the B-band ZPs so that we can correct 
each source for extinction individually. After convolving 
the images to match the largest seeing (often the V or I 
band), we applied the B- and z-band ZPs to their cor¬ 
responding photometric catalogs generated in two-image 
mode based on detections in R (see §3.1), and correct 
for extinction differentially using the method described 
above. 


3.3. Rest-frame Magnitudes and Colors 

In f[5l we use rest-frame absolute B-band magnitudes 
when selecting galaxies for the analyses, and we also use 
rest-frame U — B colors for color magnitude diagrams 
(CMDs) and B — V colors to calculate stellar masses 
in 35.31 These rest-frame magnitudes and colors are 
calculated using Q — 4 LDP reds hifts (see HH) with 
the IDL code KCORRECT (v4_2; iBlanton fc Roweisl 
I2007D . This code AT-corrects and fits the BVRIz pho¬ 
tometry (or subset thereof) using a set of templates con¬ 
structed from Bruzual-Charlot stellar evolution models 



18 19 20 21 22 23 24 25 

Rauto 


Fig. 1.— ( Top panel) Differential number counts of R-band de¬ 
tected sources per 0.1-sized magnitude bin ( dN/dm ) as as func¬ 
tion of -Rauto- We have not removed stars from the distribution. 
The solid line shows a fit to the distribution using a linear regres¬ 
sion, whileanestimate of dN/dm for the core VLT photometry 
from I White et all (120050 is shown as a dotted line. ( Bottom panel) 
Residuals from the best-fit show dN/dm follows a power law until 
Rauto = 23.3 (vertical dashed line), a clear sign of incompleteness 
beyond that magnitud e. _ 

iBruzual fc Chariot! (l2003f ). For some of our lowest red- 
shift cluster galaxies we lack observed-frame filters that 
are blueward of rest-frame U, such that KCORRECT 
must extrapolate to measure the U — B color. Based 
on their appearance in the CMDs, it is likely that their 
extrapolated rest-frame colors are too blue. By compar¬ 
ing their observed-frame V — I colors with similarly blue 
galaxies that have observed-frame B-band magnitudes, 
we find comparable V — I colors. This demonstrates that 
while KCORRECT may measure a U—B color that is too 
blue, the galaxies still meet our definition of blue galax¬ 
ies (see 35.31) . and in any event, removal of these galaxies 
with extrapolated rest-frame colors does not fundamen¬ 
tally alter any of our main conclusions. 

3.4. Photometric Completeness 

We estimate our magnitude-limited completeness by 
examining the galaxy number counts as a function of 
magnitude (Figure [T|. Differential number counts with 
magnitude (dN/dm) follow a power-law distribution un¬ 
til the shape of the curve turns over once the cata¬ 
log starts to become incomplete, with deeper catalogs 
turning over at fainter magnitudes (e.g., Figure 1 of 
[White et al.1 120051 although they applied aperture cor¬ 
rections to their magnitudes that result in a slightly dif¬ 
ferent slope and sharper cutoff at the faint end). In Fig¬ 
ure [U we fit dN/dm using a linear regression, and find it 
follows a power-law until Rauto ~ 23.3, at which point 
the distribution turns over. This turnover happens at a 
similar magnitude in all of our fields and provides an es¬ 
timate of our completeness, i.e. we are photometrically 
complete to Rauto < 23.3. 

4. LDP SPECTROSCOPIC DATA 

In this section, we present details on the LDP target 
selection, as well as the redsliift-fitting procedure and 
results. 


4.1. Target Selection 
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We utilize the LDP and the IMACS camera on the 
Magellan I Baade 6.5m telescope at Las Campanas Ob¬ 
servatory. This instrument provides spectra with a reso¬ 
lution of 7 Z = A/AA « 20-120 from red to blue wave¬ 
lengths, an improvement over the resolution achieved 
with photometric redshifts ( 1Z ~ 5). The correspond¬ 
ing redshift precisio n is als o improved, as is the overall 
accuracy. I Coil et all (1201 il l present more details about 
the prism and camera characteristics. 

Of the original 20 EDisCS fields, four were not 
targeted with the LDP, and so do not appear in 
this paper. C11119—1129 and C11238—1144 do not 
have near-infrared (NIR) data; the former contains 
a a = 166 km s -1 cluster while the latter has 

only fou r spectroscop ic reds hifts at the cluster distance 
(iMilvang-Jensen et al.1 I2008H . C11122—1136 does not 
contain a confirmed cluster, and 01202—1224 was not 
observed owing to the limited telescope time available. 

We obtained the LDP data during two observing runs, 
from 7-9 February 2008 and 27 30 March 2009 (Table 2). 
Slit dimensions are 1" x 0.8", com pared to 1" x 1 .6" for 
the bulk of the PRIMUS survey (iCoil et al.l 1201 Ilf : this 
choice allows the placing of « 1800-2800 slits per mask. 
We chose exposure times of 64 x 60 s per mask and used 
nod-and-shuffle mode to improve sky subtraction. 

We observed each held with two masks, except 01232- 
1250, which was observed with three masks. Portions of 
each Held are masked out due to the presence of bright 
stars. The FOV covers ~ 0.2 square degrees around each 
cluster, corresponding to cluster-centric distances of ~ 6- 
8 Mpc. Because each mask in a given Held has a different 
center, the final footprint for each field has a non-regular 
shape. 

There are 1-2 x 10 4 sources in our WFI catalog within 
each LDP footprint. Of these, we target « 3000-5000 ob¬ 
jects per field with the LDP (« 20% of potential targets, 
although the percentage ranges among the fields from 15 
to 40%). Galaxies are targeted depending on their 77- 
band magnitude relative to the brightest cluster galaxy 
(BCG). Priority #1 targets have Rauto.bcg — 1 < 
-Rauto ~ 23. After targeting those, additional slits (pri¬ 
ority #2) are placed on sources with 19<Rauto < 
Rauto.bcg — 1- In cases where the Rbcg.auto < 20, 
there are no priority #2 targets. Finally, we place slits 
on any galaxies that do meet these criteria but are ca¬ 
pable of being targeted in that mask; these “filler” slits 
are ~ 10% of the total. The R-band ranges of non-filler 
targets per field appear in Table 3. 

The mean separation between adjacent slits for an indi¬ 
vidual mask is ~ 20", with a minimum separation of 10". 
However, multiplexing done with multiple masks per field 
increases the sampling density, with a mean separation 
of « 15" and 15-20% of slits separated by < 10" (with 
the closest pairs « 1" apart). 

4.2. LDP Redshifts 

The PRIMUS reduction pipeline simultaneously fits 
the spectral and photometric data to a set of galaxy 
templates at different redshif ts an d calculates a best-fit 
X 2 value at each redshift (iCool et al.ll2013lh While the 
relative astronometry of our observing masks was accu¬ 
rate, they were mildly offset in absolute astronometry. 
We calculated a new astrometric solution for the WFI 


TABLE 2 

LDP Observing Log 


Run* 

Cluster 

Seeing (") 

i 

C11040.7-1155 

0.7 

i 

01054.4-1146 

2.2 

i 

01054.7-1245 

0.6-1.1 

i 

01103.7-1245 

0.6-1.0 

i 

01216.8-1201 

0.7 

i 

01227.9-1138 

0.5-0.9 

2 

01018.8-1211 

0.6-0.7 

2 

01037.9-1243 

0.5-1.0 

2 

01059.2-1253 

0.5-0.6 

2 

01138.2-1133 

0.5-0.6 

2 

01232.5-1250 

0.4-0.6 

2 

01301.7-1139 

0.5-0.7 

2 

01353.0-1137 

0.4-0.6 

2 

01354.2-1230 

0.4-0.6 

2 

01411.1-1148 

0.4-0.7 

2 

01420.3-1236 

0.5-0.7 


Run 1 took place 6 Feb 2008 to 8 Feb 
2008; Run 2 took place 27 Mar 2009 to 
30 Mar 2009. 


TABLE 3 

Photometric Targeting Criteria 


Cluster 

Rauto.bcg 

Priority #1 

Priority #2 

01018.8-1211 

19.64 

18.64-23.13 


01037.9-1243 

19.79 

18.79-22.89 


01040.7-1155 

21.17 

20.17-22.99 

18.99-20.17 

01054.4-1146 

21.20 

20.20-22.93 

18.93-20.20 

01054.7-1245 

21.09 

20.09-22.97 

18.97-20.09 

01059.2-1253 

19.20 

18.20-23.11 


01103.7-1245 

22.87 

21.87-22.99 

18.99-21.87 

01138.2-1133 

20.03 

19.03-22.67 


01216.8-1201 

20.56 

19.56-23.07 

19.07-19.56 

01227.9-1138 

21.06 

20.06-22.99 

18.99-20.06 

01232.5-1250 

19.12 

18.12-23.05 


01301.7-1139 

19.56 

18.56-23.05 


01353.0-1137 

20.29 

19.29-23.08 

19.08-19.29 

01354.2-1230 

21.27 

20.27-22.97 

18.97-20.27 

01411.1-1148 

20.79 

19.79-23.08 


01420.3-1236 

20.09 

19.09-22.95 



Note. — Priority #1 corresponds to 

^ AUTO, BCG - 1 < Aauto <23; Priority #2 corresponds 

1 9<R< -Rauto.bcg — 1, when Rbcg.auto > 19- 


imaging after the LDP slit positions were determined by 
cross-correlating the WFI catalog with the slits using a 
1" matching threshold. The pipeline treats photometric 
data similarly to one pixel of the spectrum. The chi 2 
fit is determined primarily from the LDP data, but the 
photometric data helps distinguish between redshift solu¬ 
tions. On the basis of the y 2 distribution, both a best-fit 
redshift and a redshift confidence parameter, Q are cal¬ 
culated. The integer values of Q range from 2 to 4, with 
Q = 4 being the most secure. Example spectra of four 
cluster galaxies appear in Figure [21 Furthe r deta ils on the 
redshift-fitting procedure appear in ICool et al.l (120131) . 

For objects with more than one redshift measured 
(given the multiple masks), we take the redshift with 
the higher quality flag, Q. In cases where there are mul¬ 
tiple redshifts with the same Q, we randomly select one. 
This is done to avoid averaging significantly discrepant 
redshifts when they exist (see Figure [3] for outlier rates.) 

We define cluster membership as galaxies that have a 
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Fig. 2.— Sample LDP spectra with Q = 4 for galaxies in four clusters, spanning 21.3 < Rauto < 22.6 and < 0.4 < z < 0.7. Purple 
squares and red circles are the LDP spectra values in each slit of the nod & shuffle, while the continua spectra shown are those of the 
best-fit templates. The spectra have been normalized to unity at 6800A. Typical errors on the LDP data are shown at the bottom right of 
each panel. Also plotted are magnitudes ( VRI , BVRIz, BVRIz , and BVRI , respectively), and prominent spectral features are shown as 
vertical dashed lines. The LDP redshifts agree with those measured from FORS2 to Sz = 0.0024, —0.0132, —0.0139, and 0.0012. 


Q = 4 best-fit redshift within ±0.02 (±6000 km s _1 ) of 
the cluster redshift, which is approximately three times 
the precision of the LDP redshifts (see M.2.11) . We choose 
a fixed cut in redshift, rather than a multiple of the 
cluster a, because the LDP uncertainty is larger than 
any velocity dispersion in our sample. We do not make 
any spatial cuts, because we are interested in galaxies at 
large cluster-centric radii. This selection results in 1,602 
galaxies that we place in the cluster environment. 

A summary of the number of LDP targets, redshifts, 
and cluster members is presented in Table 4. Two of 
the LDP-observed clusters do not appear in this study. 
C11103.7-1245, at z = 0.95, with only 3 LDP-selected 
cluster members, does not have enough cluster members 
for a meaningful analysis. 01138.2-1133a, at z = 0.4548, 
lies too close to the redshift of 01138.2-1133, z = 0.4796, 
to distinguish between members using the z c ius±0.02 se¬ 
lection. We therefore include the latter cluster, with the 
caveat that some contamination may come from galaxies 
belonging to the former. We note that 01138.2-1133 is 
not a significant outlier in any of the analyses that follow. 

We include a catalog with this paper containing the 
cross-matched WFI BVRIz photometric sources, LDP 
spectroscopic redshifts, and VLT/FORS2 spectroscopic 
redshifts, as well as rest-frame UBVRIz magnitudes and 
stellar masses (see §5.3). 

4.2.1. LDP Redshift Precision 

We assess the precision of our LDP-derived redshifts 
(~ldp) by comparing them to the subset of 427 galax¬ 


ies also observed with VLT/FORS2 (^spec) over a 
wide range of redshifts and with p hotom etric redshifts 
(-zphot) calculated in iPello et ahl (l2009f l from BVIK, 
BVIJK , and VRIJK imaging of the cluster cores; these 
filter combinations were chosen based on the initial red¬ 
shift estimate of the cluster. We match galaxies within 
1", and show the results of these comparisons in Figure[3] 
We only consider galaxies with Zldp < 0.85, which is 
just above our highest redshift cluster (C11216.8 — 1201 
at z = 0.79); considering the full range of redshifts that 
PRIMUS fits (out to z = 1.2) affects neither the precision 
nor outlier rate significantly. 

The LDP-derived redshifts are more precise than the 
photometric redshifts by an order-of-magnitude, their 
RMS being <t(|;sldp — £spec|) = 0.006 compared to 0.08 
for the photometric redshifts. The outlier rates of LDP- 
derived redshifts, defined as |zldp — £spec| > 0.02, de¬ 
pend on the quality cut, and range from 25% (Q > 2) to 
20% (Q > 3) to 17% (Q = 4). 

In Figure [4] we plot the LDP redshift precision and 
outlier rate dependence on R-band magnitude. The pre¬ 
cision is constant, even at faint magnitudes. However, 
the outlier rate increases significantly for I?auto > 21. 
This increase reflects the larger number of Q = 2 and 
3 spectra at these magnitudes (see M.2.2\ ) . The outlier 
rate is constant with magnitude for spectra of a given 
Q value (although the number of outliers changes). We 
do not find any significant dependence of the redshift 
precision on V — I color (Figure [5]). We have the best 
statistics for the Q = 4 redshifts (which make up most of 
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Fig. 3. — (Left panels) Comparison of LDP-derived (zldp) and photometric redshifts (^phot) with spectroscopic redshifts (^spec)- The 
LDP redshifts are split showing different cuts in the quality flag, Q. The outliers in ^ldp systematically underestimate the redshift. In 
the bottom panel, open circles show all galaxies with zpnoTt while filled circles mark the subset that also have LDP redshifts (Q > 2). 
{Right panels) Histograms of the residuals from the left panels. The dotted line in the bottom panel include all galaxies with zpHOT 
measurements, while the solid line shows .zphot residuals for the subset that also have LDP redshifts (Q > 2). Vertical lines show ±0.02, 
which is the size of the redshift interval used in selecting cluster galaxies. Percentages in the upper left corner show the fraction of outliers 
outside this interval. The precision with the zldp = 0.006) is an order-of-magnitude improvement over the photometric redshifts. 


our redshifts), and find a stable outlier fraction of 17%. 

From Figures [3] and 01 it is apparent that the outliers 
in LDP redshift are skewed towards lower redshift val¬ 
ues (i.e., the LDP fits a lower redshift than the “true” 
one). Although the outlier rate does not depend on 
galaxy color, we find a significant dependence on [Oil] 
emission. Being relatively close to the Balmer break, the 
blending of the two spectral features “drags” the break 
to a lower redshift, consistent with the bias evident in 
Figure [3j When we consider Q = 4 LDP spectra, we 
find that galaxies with [Oil] equivalent widths (EWs) 
of ~ 5A have an outlier rate of 35%, more than dou¬ 
ble the outlier rate for all galaxies (17%). Moreover, at 
both higher and lower EWs, the outlier rate drops. This 
can be understood as galaxies with weaker [Oil] emis¬ 
sion not suffering from this blend “dragging” the Balmer 
break to a lower redshift, while galaxies with stronger 


[Oil] emission having fines that become the dominant 
redshift feature. Therefore, we are more likely to miss 
cluster galaxies that have modest [Oil] emission at an 
outlier rate twice as high as for the full galaxy popula¬ 
tion, or ~ 4% of the cluster sample, and have an some 
enhanced contamination from field galaxies with modest 
[Oil] emission at higher redshift. 

4.2.2. Spectroscopic Completeness 

Figure [5] shows the distribution of Rauto (analogous 
to Figure |T]) for the LDP targets with successfully mea¬ 
sured redshifts. The different curves show the distribu¬ 
tions for Q = 2, 3, and 4 redshifts. At the brightest 
magnitudes, the vast majority of redshifts have a secure 
Q = 4 flag. However, at Rauto ~ 20, redshifts with lower 
Q flags begin to appear in significant numbers. For all Q 
values, the distribution turns over before Rauto ~ 23; 
our photometric catalog is therefore complete to fainter 
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TABLE 4 
LDP Information 


Cluster 

(i) 

-^phot 

(2) 

-Ntargets 

(3) 

Wldp 

(4) 

N Q =4 

(5) 

•^memb 

(6) 

C11018.8-1211 

11868 

1454 

1245 

734 

73 

C11037.9-1243 

9569 

1971 

1781 

1050 

44 

C11037.9-1243a 





171 

01040.7-1155 

8581 

2373 

1848 

827 

29 

01054.4-1146 

11551 

2518 

2080 

1120 

43 

01054.7-1245 

8923 

2327 

1945 

995 

72 

01059.2-1253 

7553 

1998 

1792 

1247 

131 

01103.7-1245* 

8203 

2230 

1245 

522 

2 

01103.7-1245a 





18 

01103.7-1245L 





18 

01138.2-1133 

10293 

1671 

1511 

1202 

91 

01138.2-1133a* 






01216.8-1201 

10741 

2491 

2018 

972 

46 

01227.9-1138 

8710 

2564 

1898 

1137 

75 

01227.9-1138a 





93 

01232.5-1250 

9353 

2320 

2078 

1414 

158 

01301.7-1139 

10147 

1526 

1318 

968 

126 

01301.7-1139a 





136 

01353.0-1137 

13751 

1671 

1450 

1019 

28 

01354.2-1230 

12438 

2358 

2135 

1514 

28 

01354.2-1230a 





69 

01411.1-1148 

9040 

1554 

1379 

1029 

67 

01420.3-1236 

9599 

1766 

1555 

953 

84 

Total 

160320 

32792 

27278 

15828 

1602 


Note. — All numbers only include galaxies brighter than R < 
22.6, our spectroscopic completeness limit. Numbers for columns 
2-5 for serendipitously discovered clusters are suppressed as they 
are in the same field as the primary cluster. (1) cluster name; 
(2) number of photometric sources; (3) number of LDP targets; 
(4) number of successfully extracted LDP spectra; (5) number of 
Q = 4 LDP spectra ; (6) number of cluster members (defined by 
z c \ us ± 0.02). 

* Doe s not appear in the analyses presented in this paper (see 

§E3. 



Rauto 


Fig. 4.— (Central panel) Residuals between LDP and FORS2 
redshifts as a function of R-band magnitude for Q = 4 (black cir¬ 
cles), Q = 3 (blue triangles), and Q = 2 (red squares) redshifts. 
(Top panel) Histograms of R-band magnitude and zpj^p — .zspec 
residuals. (Bottom panel) Outlier rate (|A 2 :| > 0.02) as a func¬ 
tion of R-band magnitude. The outlier rate is approximately flat 
for a given Q flag. The vertical dotted line is our spectroscopic 
completeness limit. 



0.0 0.5 1.0 1.5 2.0 2.5 

(V-I) 


Fig. 5.— ( Central panel) Residuals between LDP and FORS2 
redshifts as a function of V — I color. Symbols are the same as 
Figure |4] (Top panel) Histogram of V — I color. (Bottom panel) 
Outlier fraction as a function of V — I color. While there is a 
bump in the outlier fraction of Q = 4 redshifts at (V — I) ~ 1.2, 
no significant trends with V — I across the full color range exist. 

magnitudes than our spectroscopic one. The full distri¬ 
bution (including all Q values) departs from a power-law 
at i?AUTO ~ 22.6, which we take as the estimate of our 
spectroscopic completeness. 

4.2.3. Radial Completeness 

We also quantify the percentage of sucessfully mea¬ 
sured redshifts as a function of clustercentric distance. 
We consider the fraction of sucessfully measured red¬ 
shifts relative to the number of photometric sources, re¬ 
stricting both to galaxies brighter than our spectroscopic 
completeness (i?AUTO < 22.6), as a function of angular 
distance from the cluster (d c ius)- Figure [7] shows that the 
percentage of targets with measured redshifts is ~ 20- 
30%, depending on Q-cut, out to ~ 10'. Converting this 
to a physical distance for our typical clusters puts the 
drop-off at ~ 4 Mpc. The percentage then drops off, as 
most fields have only one mask coverage at these radii, 
and approaches zero smoothly rather than abruptly be¬ 
cause the offset placement of masks leads to an edge that 
is not spherically symmetric about the cluster. 

Given the surface density of slits on the sky, we are also 
relatively insensitive to close pairs. While multiplexing 
with 2-3 masks per field allows us to measure redshifts 
for galaxies ~ 1" apart that are not near the edges of the 
footprint, this distance is small compared to the average 
separation between adjacent slits (« 20"). Only « 10% 
of slits have separations of 10" or less. Compared to fre¬ 
quency of slits with separations of 10-30", the frequency 
of slits separations < 10" is only ss 40% of that value. 
This affects our ability to find cluster galaxies within the 
cluster core, where the galaxy surface densities are higher 
(see (J5|). 

5. CLUSTER INFALL REGIONS 

We use the the ory of secondary infall 
dFillmore fc Coldreichl 11984 iBertschinged 1198.4 
iWhite fe Zaritskvl Il992f) to estimate (1) the masses 
that our clusters are expected to accrete by z = 0 
(Afinfaii) and (2) the projected radii at the cluster 
redshifts that encloses Mi n f a n (i?i n f a ii). Using the latter 
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Fig. 6. — (Top left panel) Histogram of i^-band magnitudes for 
LDP-targeted galaxies with bin sizes of 0.1 for different values of Q. 
The spectroscopic completeness limit at Rauto — 22.6 is shown 
as a vertical line. (Bottom left panel) The fraction of targets with a 
successfully extracted spectrum, i.e. number of spectra divided by 
number of slits, as a function of magnitude. At bright magnitudes, 
most of the successfully extracted spectra have Q = 4. (Right 

panels) Same as left panels, only as a function of (V — I) color. 



Fig. 7.— Fraction of photometric sources with LDP redshifts 
(^ldp/^phot) as a function of cluster-centric angular distance 
(d c ius)- A distance of 10 ; corresponds to ~ 4 Mpc, depending on 
the cluster redshift. Bin sizes are 0.5 arcminutes, and solid, dotted, 
and dashed lines include redshifts with Q = 4, Q > 3, and Q > 2, 
respectively. 

with the LDP data, we calculate the number of galaxies 
in the infall region (N in { aU ). In 115.ll we calculate the 
expected evolution of the clusters in terms of mass and 
compare to semi-analytic models. In H5.2I we focus on 
the dependence of Mnfaii with clu ster velocity dispersion 
and quantify its scatter. In 115.31 we examine the 
quiescent fraction of galaxies in different environments, 
using the red-sequence galaxy fraction as a proxy, and 
quantify the amount of clustering in the infall regions. 

5.1. The Secondary Infall Model 

While previous studies of cluster infall regions have 
used the caustic technique (e.g., iGeller et all 119991 
Rines et all 120031 : iRines fc Diaferiol l2006t IS err a et ahl 
201 111 , which identifies curves in galaxy position-radial 
velocity phase space that encompass those galaxies that 
are gravitationally bound to the cluster, such an analysis 


requires redshift precision greater than that of the LDP. 
Alternatively, we identify the infall regions of our clusters 
using the theory of secondary infall to define projected 
radii that encompass the infall region. 

The secondary infall model describes how shells of mass 
centered on a cosmic perturbation evolve over time. The 
shells begin by expanding outward, until a time R urn 
when they turn around due to the pull of gravity. The 
shells do not cross during this time, and there is a crit¬ 
ical mass, M*, enclosed by the shell that is marginally 
bound. All shells enclosing a mass less than M* even¬ 
tually turn around at different times and collapse, while 
shells at larger radii continue to expand foreve r. We fol¬ 
low the equations of lWhite fe Zaritskvl (11992D . who as¬ 
sume an open universe with Da = 0. Keeping Da = 0 
when calculating M; n f a n and Rinfaii does not significantly 
affect the an alysis given the physical scales involved (e.g., 
iDel Popolo et al.ll2012j ). as we also confirm with our own 
comparison to cosmological simulations discussed below. 
For determining global quantities (e.g., connecting a time 
with a redshift and age), we continue to use Da = 0-7. 

Beg inning with Equation 13 of iWhite fc Zaritskvl 

(EMI, 


t 


turn 


(Me nc) 


7r D 

2 H(1 — D) 3 / 2 



-3/2 


(3) 

we solve for the critical mass M* of each cluster by setting 
the turnaround time Rum to half the age of the Universe 
at z = 2 i c i us and the enclosed mass (M enc ) to M 200 (from 
Table 1). For this value of M*, the equations of secondary 
infall pre dict the virial mass co rres ponding to the a mea¬ 
sured in lHallidav et all (T2004 ) or IMilvang-Jensen et all 
<2QQ3). Once M* is determined, we use Equation 2 with 
Rum equal to half the present age of the Universe, to de¬ 
termine the predicted M enc at 2 : = 0 (M 2 oo,z=o); a shell 
turning around at this time is the outermost mass shell 
that will reach R = 0 at 2 = 0. The turnaround radius 
(Rturn) for this shell is related to R ur n through the simple 
equation for a free-falling test particle, 


Rurn(ARoO,z=o) 


/? 3 

x Hurn 


2GM- 


200,z=0 


(4) 


While we solve for iRurn for each cluster, which turns 
around when the Universe is half its present age, what 
we are truly interested in is the location of these shells 
at 2 = 2 c ius- Therefore, we use the equation of motion 
for a uniform mass shell to evolve iRurn to 2 = 2 c i us . 
These radial distances are the infall radii (iRnfaii), the 
outer boundaries of the relevant infall regions for 2 = 0 
observations. 

In defining radial distances, we center on the location of 
the BCG. However, the BCG may be offset from the dis¬ 
tribution of mass, which would affect the definition of the 
infall region. We estimate the magnit ude of these offsets 
from Figure 6 of 1 White et al.1 (12005( 1. which marks the 
BCG position relative to adaptively smoothed contours 
of cluster-galaxy surface density. The offsets are < 10% 
of Rinfaii for all clusters except C11037-1243, whose BCG 
is offset by ~ 25% of Rafail- These v alues are larger tha n 
the typical offsets found at 2 ~ 0.5 bv lZitrin et all (l2012f) . 
but they find the offsets are positively correlated with 
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redshift and our clusters lie at higher z than the ir sam ple. 
In addition, the galaxy distributions from [White et al.1 
(120051) include galaxies with photometric redshifts con¬ 
sistent with being close to the cluster redshift, and there¬ 
fore include a non-negligible number of interlopers that 
make the centering less precise. If we define the cen¬ 
ter to be the mean right ascension and declination of 
the VLT/FORS2 spectroscopic sources, then 41 galaxies 
(12% of the infalling population) are either removed or 
added by the new definition. However, the total number 
of infalling galaxies changes by less than 2% (because 
some are added, while some are removed from the infall 
region), and the fraction of red galaxies 1( 15.31) changes by 
only 0.4%. As a final check, we randomly apply offsets of 
~ 0.2 Mpc in various directions from the BCGs and find 
similar numbers of galaxies are affected by the redefined 
infall regions. We conclude that reasonable uncertainties 
in the centering of the clusters do not strongly impact 
our results or conclusions. 

In Table 5, we present the results of the models for 
our clusters. In addition to i?i n faii, we calculate the pre¬ 
dicted mass (and corresponding velocity dispersion) at 
z = 0 for our sample. The infall radii range from 1.2 
to 6.7 Mpc; the ratio of AWaii to i ?200 is set entirely by 
the redshift of the cluster, in that higher -z clusters have 
larger ratios, and range from 3.0 to 4.2i? 2 oo- These are 
smaller than the “turnaround radii” calculated in other 
studies of cluster infall regions, such as those found in 
dRines fe Diaferiol2006l ~ 4.75i? 2 oo)- However, their def¬ 
inition of infall region includes all galaxies that, with a 
velocity less than the cluster escape velocity, will eventu¬ 
ally become incorporated into the cluster given enough 
time, while our model only includes galaxies that could 
have reached the center of the cluster by z = 0. 

We compare our i nfall radii to models using the Mil¬ 
lennium Simulation (Springel et al.ll2005h , estimating the 
infall radii in the simulation using the fraction of galax¬ 
ies at a given clustercentric distance that come to lie 
within R 200 at z = 0. We considered 174 ~ 1O 14 M 0 
halos at z ~ 0.6, which have i?i n faii ~ 3.5i?200 according 
to our analytic modeling. We find that 67 ± 7% of the 
galaxies within 3.5i?200 at 2 = 0.62 ultimately lie within 
the virial radius of the descendant halo (the errors are 
the 15-85 percentiles). Because we assume that 100% of 
galaxies within i?i n faii become “cluster galaxies” of the 
descendant halo, our prediction is good to « 33%. The 
difference arises from a variety of effects, including the 
assumption of spherical symmetry, the definition of an 
infalling galaxy as one that reaches R = 0 at z = 0, 
and the presence of “backsplash” galaxies, which pass 
through the clust er core an d then continue out to radii 
larger than i? 2 oo dCfill et, al.ll2005 ). 

Our models predict an increase in cluster mass of 26- 
57% from the observed epoch to the current one (col¬ 
umn 7 of Table 5). Velocity dispersion increases as 

<7ocM 1 / 3 [n A + n 0 (l + 2) 3 ] 1/6 , (5) 

which means the mean growth corresponds to a increas¬ 
ing by ~ 15-25%. The predicted a at z = 0 agrees to 
within ~ 1 0% o f the predictions at a given mass that 
iPoggianti et al.l (|2006t l computed for 2 = 0.6 clusters 
by com bining the high-resolution N-body simulations 
of IWechsler et al.l (|2002i l with cluster concentration 


parameters from iBullock et all (120011) . 


5.2. Number of Infalling Galaxies 

We estimate the richness of the infall regions (R 200 < 
R < i?i n faii), which we define as the number of clus¬ 
ter galaxies in the infall region (W n f a n) above an ab¬ 
solute R-band magnitude of Mb = —19.2 (correspond¬ 
ing to Rauto ~ 22.6 for our highest redshift cluster at 
z = 0.79). 

Because we do not have redshifts for every galaxy 
above this magnitude limit, to estimate -/Vi n f a ii we use 
the number of photometric sources in the radial range 
R 200 < R < Rinfaib Al p hot, multiplied by an estimate 
of what fraction of N p hot lies within A z ± 0.02 of the 
cluster. We determine this fraction from the ratio of 
LDP-selected cluster members, M mem b, to the number of 
LDP slits, ATgiits, limiting both to the infall region. This 
procedure accounts for the incomplete spatial sampling 
due to chip gaps and masked bright stars. We estimate 
the contamination from field galaxies (/contain) by us¬ 
ing the fraction of field galaxies at that cluster redshift; 
these galaxies he at z c i us but are observed in fields other 
than that particular cluster (excluding any that lie within 
another EDisCS cluster; see d5.3l) . Therefore, AVifaii is 
calculated as 


Afinfall = A^phot ^ emb (1 — /contain)- 


slits 


( 6 ) 


We also estimate the number of cluster galaxies within 
i ?200 at z = Zcius, Vciuster, using the same methodology, 
but apply an additional correction to account for close 
pairs that the LDP might miss (cf. ( 14.2.31) . The fre¬ 
quency of LDP slits with separations between 10-30" is 
fairly constant, changing by less than 15% over this sep¬ 
aration range, but below ~ 10" it plummets to « 30% of 
that frequency. Given that 35% of cluster members iden¬ 
tified with VLT/FORS2 have such separations or less, 
one way to correct for this would be to divide N c i usteT by 
the reduced frequency of sampling (30%) as well as one 
minus the fraction of close pairs (35%); the correction 
would therefore involve dividing by 20%. Alternatively, 
we can look at how many of the core cluster galaxies 
targeted with VLT/FORS2 were also targeted with the 
LDP and divide by that fraction, which is 18%. Because 
the correction factors are comparable and the latter is 
derived more directly, we use 18% as the corrective fac¬ 
tor. Although applying this correction will add uncer¬ 
tainty to the normalization of N c i us t el7 our derived values 
of V c i us ter agree within the uncertainties to previous es¬ 
timates made using spectrosco pically-confirmed EDisCS 
galaxies IPoggiant i et al.ll2010H . 

In Figure 0 we present a plot of V c i uste r and A%f a n 


versus a. Errors on V c i uste r and A r ; n f 
nian (calculated using the equations of 

a .n are Poisso- 

Gehrcls 

1983) 

and errors in a come from Hallidav et al. (2004 

) and 


iMilvang-Jensen et al.1 (|2008t ). We perform a linear re¬ 
gression that follows a Bayesian approach and accounts 
for errors in both a and N , usin g the IDL routine 
LINMIX_err.pro written by iKellvl (l2007t ). We find 
Cluster OC (j 21±10 and A/infaii oc (T 1 - 6 ^ 1 - 0 . While the 
trend for V c i uste r has a steeper depen dence on a tha n that 
of the z ~ 0 clusters observed bv lFinn et al.1 (|2008l h who 
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TABLE 5 

Mass Infall Model Results 


Field 

(1) 

Cluster 

(2) 

z 

(3) 

-^infall 

(4) 

■^infall 

?sT 

Af200,z=0 

(6) 

M200, z =0 

M^o 

& z=0 
(8) 

o~z=0 

cr 

(9) 

1 

01018.8-1211 

0.4734 

3.06 

3.28 

2.01(14) 

1.32 

533 

1.10 

2 

01037.9-1243 

0.5783 

2.07 

3.60 

5.67(13) 

1.40 

357 

1.12 

3 

C11037.9-1243a 

0.4252 

3.31 

3.12 

2.72(14) 

1.28 

584 

1.09 

4 

01040.7-1155 

0.7043 

2.75 

3.94 

1.27(14) 

1.50 

478 

1.14 

5 

01054.4-1146 

0.6972 

3.87 

3.92 

3.55(14) 

1.49 

673 

1.14 

6 

01054.7-1245 

0.7498 

3.32 

4.05 

2.21(14) 

1.53 

581 

1.15 

7 

01059.2-1253 

0.4564 

3.19 

3.23 

2.33(14) 

1.31 

558 

1.09 

8 

01103.7-1245a 

0.6261 

2.20 

3.73 

6.61(13) 

1.43 

379 

1.13 

9 

01103.7-1245b 

0.7031 

1.66 

3.93 

2.78(13) 

1.49 

288 

1.14 

10 

01138.2-1133 

0.4796 

4.62 

3.30 

6.88(14) 

1.32 

804 

1.10 

11 

01216.8-1201 

0.7943 

6.70 

4.15 

1.82(15) 

1.57 

1183 

1.16 

12 

01227.9-1138 

0.6357 

3.76 

3.76 

3.29(14) 

1.44 

648 

1.13 

13 

C11227.9-1138a 

0.5826 

2.21 

3.61 

6.93(13) 

1.40 

382 

1.12 

14 

01232.5-1250 

0.5414 

6.95 

3.49 

2.21(15) 

1.37 

1199 

1.11 

15 

01301.7-1139 

0.4828 

4.34 

3.31 

5.69(14) 

1.33 

755 

1.10 

16 

01301.7-1139a 

0.3969 

2.37 

3.02 

1.05(14) 

1.26 

423 

1.08 

17 

01353.0-1137 

0.5882 

4.33 

3.63 

5.16(14) 

1.41 

746 

1.12 

18 

01354.2-1230 

0.7620 

4.27 

4.08 

4.70(14) 

1.54 

749 

1.16 

19 

01354.2-1230a 

0.5952 

2.82 

3.65 

1.42(14) 

1.41 

486 

1.12 

20 

01411.1-1148 

0.5195 

4.54 

3.43 

6.27(14) 

1.35 

785 

1.11 

21 

01420.3-1236 

0.4962 

1.38 

3.35 

1.82(13) 

1.34 

240 

1.10 

Note. 

— (1) cluster field; 

(2) cluster 

name; 

(3) cluster redshift; i 

(4,5) infall radius in 

units of 


Mpc and units of observed-epoch virial radii; (6,7) virial mass evolved to z — 0 in units of M©(with 
powers of ten in parentheses) and units of observed-epoch virial masses; (8,9) velocity dispersion 
evolved to z — 0 in units of km s —1 and units of observed-epoch a. 


found iV c i us ter oc <7 1-4 , the results are within the uncer¬ 
tainties. Given the large errors, we find no statistically 
significant difference in the intercepts in either region. 
While more precise measurements are needed to reduce 
the uncertainty, our best-fit scalings between iV c i U ster and 
Ajnfaii with a are comparable. 

Given the similar slopes, we now examine whether clus¬ 
ters of greater mass accrete proportionally more or fewer 
galaxies over this redshift interval. Such behavior would 
have ramifications for the cr-d ependent increase in SOs as 
a fraction of cluster galaxies (jJust et al.ii2010i ). For ex¬ 
ample, if more-massive systems accrete a larger percent¬ 
age of their galaxies at late times relative to less-massive 
systems, then it could be that the (proportionally larger) 
infalling population diluted any increase in the SO frac¬ 
tion in these systems, rather than that the less-massive 
systems are intrinsically more efficient at converting spi¬ 
rals to SOs. In Figure [51 we compare the ratio of 7Vi n f a n 
to lV c i us ter as a function of cr. The effects of merging be¬ 
tween galaxies, which would decrease this ratio assuming 
no mergers occur within the virial radius, where rela¬ 
tively velocities between galaxies is much higher, is also 
shown. We find that the relative size of the infalling pop¬ 
ulation does not scale with cluster velocity dispersion. 
However, there is a considerable range in the ratio, from 
« 15-300%, with typical values between ~ 30-200%. Al¬ 
though most of this is due to the scatter in -/V c i uste r (see 
Fig [8j), this result highlights the importance of account¬ 
ing for the large variation in cluster properties. 

We now use -/V c i us ter and iVi n f a n to predict the mass 
evolution of our clusters based on the LDP data. We 
account for the enhanced clustering of galaxies relative to 
the underlying mass distribution, parameterized by erg, 
the RMS fluctuation of galaxies in an 8 hr 1 Mpc sphere 
relative to fluctuation in mass, and adopt a recent value, 



Fig. 8.— {Top panel) Number of galaxies within the virial radius 
(iVciuster) as a function of cluster velocity dispersion (cr). The 
dashed line is the best-fit to the data. (Bottom panel) Same as top 
panel, only for galaxies in the infall region (N[ n f a n). 

crg = 0.81 (j.Tarosik et al.ll201lb . 

From Table 5, we find that the typical increase in mass 
predicted by our adopted secondary infall models is 26- 
57%, or in terms of galaxy number, 33-74%. This is 
shown as a gray band in Figure [9l where it is consis¬ 
tent with our measured values of -/Vi n f a n/-/V c i us t e r that are 
typically ~ 30-150% (24-110% in mass). These values 
are lower than the factor of two mass increa s e pred icted 
for 0.1 < z < 0.3 clusters of iRines et all (|2013T ). al¬ 
though their use of the caustic technique means their 
prediction should be higher that ours, because theirs 
is for the final cluster mass in the distant future, not 
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Fig. 9. — Ratio of infalling galaxies to cluster galaxies as a function of velocity dispersion (<r). The horizontal gray bar comes from 
converting the expected mass increase based on the secondary infall models (Table 5) to galaxy number using as (see text). No trend is 
present, such that clusters over this range of a accrete proportionally similar numbers of galaxies as they evolve to z = 0 (although they 
typically range from ~ 0.3-2.0). Vertical dashed lines show (from left to right) how the ratio changes if 7Vi n f a n is reduced by 20%, 40%, 
60%, and 80%, respectively, due to mergers. 


z = 0. Similarly, [Dressier et al.1 (|2013i l find that among 
very rich clusters at z ~ 0.4, the number of galaxies in 
infalling groups will roughly double the mass of the clus¬ 
ters by the present, which is larger than our estimate 
but still consistent within the scatter. That we find a 
comparable mass increase among some of our clusters 
that have lower velocity dispersions than their sample 
(« 800-1100 km s” 1 ), further supports our conclusion 
from Figure [9] that more massive systems do not accrete 
proportionally more galaxies as they evolve. 

5.3. Optical Properties of the Cluster and Infalling 
Galaxies 

To estimate the quiescent fraction of galaxies in dif¬ 
ferent environments, we use the fraction of optically-red 
galaxies. We examine the CMDs of our clusters and com¬ 
pare the red fractions of core, infalling, and field galax¬ 
ies. To construct field samples for each cluster, we select 
galaxies at the same redshift but observed in fields other 
than the cluster’s (excluding any that overlap in redshift 
with the EDisCS cluster of that particular field). Note 
that (1) this means there are significantly more galaxies 
in a given field sample than in the corresponding clus¬ 
ter, because they are drawn from multiple fields, and (2) 
we combine a subset of these field samples in some of 
the analyses below, so we distinguish between “individ¬ 
ual field samples” and a “combined field sample”; the 
reason and details for which are described below. 

For our CMDs, we use rest-frame U — B colors and 
absolute B-band magnitudes calculated using the IDL 
code KCORRECT (see ^3.31 iBlanton fe Roweii 120071) . 
In FigurcflOlwe present observed-frame CMDs for the 21 
clusters. We find well-defined color-magnitude relations 


(CMRs), where old, passive galaxies at that redshift lie, 
for all of our clusters. We measure the CMRs by assum¬ 
ing zero slope and fitting the WFI U — B colors of the 
subset of cluster galaxies that have FORS2 spectra show¬ 
ing no [Oil] emission. The CMRs measured this way are 
in agreement with the apparent red sequences of LDP- 
selected cluster galaxies (Figure [Toll. In what follows, we 
define galaxies with colors within 0.2 of the CMR or red¬ 
der as red , while the remaining galaxies are classified as 
blue. To complement the CMDs, we plot spatial maps 
of our clusters in Figure [H] with red and blue galaxies 
highlighted and the virial and infall radii overplotted. 

Because environmentally-driven galaxy evolu¬ 
tion is correlated with the_velo city dis p ersion 

of the group/cluste r _ (e g . iPo ggianti et, all 120091 : 

IBlanton fc Moustakasl 12009 : I.Tust et all 120101 )in Fig¬ 
ure [12] we plot the red fraction of the cluster and infalling 
samples for each cluster as a function of cr. To compare 
to the field, we create a “combined field sample”. The 
“individual field samples” in Figure [TO] have galaxies 
common to more than one sample. To avoid this multi¬ 
ple counting in the “combined field sample”, we select 
the “individual field samples” of 9 clusters that span the 
full redshift range from 0.4 to 0.8 but do not overlap in 
redshift. The red fraction for the combined field sample 
is also shown in Figure [12] We do not find a correlation 
between red fraction and a i n the cluster environment , 
consistent with the findings of lValentinuzzi et al.l ( 2011 s ), 
and find no correlation in the infalling sample, either. 
However, the red fraction increases as one moves from 
the most isolated environment to the cores of clusters, 
from 37 ± 1% in the field, to 43 ± 2% in the infall region 
and 63 ± 4% in the virial region. 
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Fig. 10. — Color-magnitude diagrams (CMDs) for our clusters in rest-frame U — B vs. absolute 5-band magnitude. The left frames show 
LDP-selected cluster members as circles that are colored based on our red/blue definition, while the right frames show the individual field 
samples. Dotted lines mark the cluster color-magnitude relations (CMRs). Note that because of the way field samples are constructed, 
they contain many more galaxies than their corresponding cluster. 



-22 -20 -22 -20 -22 -20 -22 -20 -22 -20 -22 -20 


Mb 


Fig. 10. — Continued. 
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Fig. 11. — Spatial maps of our clusters. Red and blue circles mark cluster members meeting our red/blue definitions, and black circles 
corresponding to the virial and infall radii are also shown. 
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Fig. 11.— Continued. 

The red fraction in the infall regions is elevated rela¬ 
tive to the field at 2.8er (98.5%) significance, although not 
as high as the fraction within the virial radius, which is 
4.7cr higher than the infall region and 6.7<r higher than 
the field. This suggests the quenching of star forma¬ 
tion begins outside the virial radius, although some of 
these galaxies may have alr eady passe d through the clus¬ 
ter virial region in the past dGill et al.ll2005l see below for 
our treatment of “backsplash” galaxies). Because galax¬ 
ies move ~ Mpc distances over ~ Gyr timescales, the 
quenching of star for matio n coul d be gin to occur in the 
infall region fe.g.. iBalogh et al]l2000fl . or even primarily 
occur there, with the higher red fraction within i ?200 due 
to the lag between the start of quenching and the time 
for its effects to become apparent. With high-resolution 
imaging, one would be able to assess whether the infalling 
red galaxies exhibit early-type morphologies or perhaps 
a transitory phase as p assive disks. 

IRudnick et al.1 (1200911 find that the total light on the 
red sequence for 16 of the EDisCS clusters must increase 
by a factor of ~ 1-3 by 2 = 0. We predict that the 
clusters in our sample will grow by a factor of ~ 1.8 
in number of galaxies (Figure [9]). Given that the red 


fractions within the cluster and infall regions are 63% 
and 43%, respectively, passive galaxies already identified 
as such in the infall regions will increase the z = 0 total 
red sequence light by a factor of ~ 1.6. While t owards the 
low end of the range predicted bv IRudnick et all ( 2009 1. 
their study include galaxies fainter than in our analysis. 
We therefore conclude that significant further quenching 
of blue galaxies in the infall regions as the clusters evolve 
to z = 0 is not required for our sample. 

“Preprocessing” has been suggested as a way of 
transforming galaxies in locally overdense clumps 


prior to their incorporation into the cluster 

(e.g., 

Zabludoff & Mulchaevi 1998!; Moran et al. 

2001$ 

Kautsch et al. 20081: Dressier et al. 2013: Haines et al.1 

2013): Lopes et al.1 2013}: Gvbulski et al.1 2014(1. 

The 


elevated red fraction in the infall region, where locally 
overdense clumps are expected to exist, is consistent 
with preprocessing. However, this result does not rule 
out an additional global mechanism for quenching 
star formation, one that affects all galaxies at a given 
clustercentric radius equally. To explore this scenario 
further, we measure the amount of clustering among 
the galaxies in the infall regions, which will provide 
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Fig. 12.— Red fraction as a function of cluster velocity dispersion for the cluster and infalling populations. Individual clusters are 
shown as open circles (with error bars suppressed for clarity), while large filled circles are bins consisting of clusters with a < 500 km s —1 , 
500 < a < 1000 km s~ 1 , and a > 1000 km s —1 . Horizontal dotted lines show the overall red fraction ±1 standard deviation. Similarly, the 
overall red fraction of the combined field sample is shown as a shaded horizontal line, with a width corresponding to ±1 standard deviation. 
There are no statistically significant trends of red fraction with a in either environment. However, the red fraction decreases as one moves 
from the cores of clusters to the field. 


more direct evidence for the association of preprocessing 
with local overdensities. We measure the fraction of 
infalling red/blue galaxies with at least one infalling 
neighbor of similar color within a projected distance 
dg, which we denote F(< dg), for values of dg ranging 
from 0.15 to 0.35 Mpc. A control sample is constructed 
similar to the “combined field sample” but restricted to 
cluster-centric radii that correspond to the infall region 
of that particular field. 

We present these distributions in Figure [13] In both 
the infalling and control samples, a higher fraction of red 
galaxies have red neighbors than blue galaxies with blue 
neighbors. The latter value in both samples agree within 
the error bars, leading us to conclude that the infalling 
blue galaxies are not significantly more clustered than the 
control blue galaxies. However, the fraction of infalling 
red galaxies with close neighbors is consistently higher 
than that in the red control sample by A F(< dg) ~ 3- 
8%, although only at 1.5-2.3cr significance. The fraction 
of infalling galaxies of either color with a close neighbor 
of either color is also higher in the infall regions than 
in the control sample by AF(< dg) ~ 3-7% (at 1.6- 
2.7a). Alternatively, we measure the red fraction among 
close pairs of galaxies and compare to those without a 
neighbor (Figure[lT]). Clustered galaxies are significantly 
more likely to be red than those without a neighbor, and 
this effect is more significant in the infall regions than 
in the control sample. Evidently, at these length scales 
the infall regions show signs of enhanced clustering of 
red galaxies, consistent with “preprocessing”, in which 
local overdensities quench star formation prior to their 
incorporation into the cluster. 

We see the enhanced clustering and elevated red frac¬ 
tion among the infalling sample ('Figures [T51 and [111) rela¬ 
tive to the control (i.e. field) sample. But despite our def¬ 
inition of clustering, based on having a neighbor within 
dg, being the same for both samples, the red fractions 


are different. Therefore, either the infalling galaxies lie 
in higher local overdensities, or they experience an ad¬ 
ditional affect unrelated to local density. The spatial 
sampling rate of our LDP spectra means that we can¬ 
not directly compare the numbers of spectroscopically- 
confirmed galaxies within dg between the two samples. 
Instead, we compare (1) the number of photometric 
sources and (2) the total i?-band luminosity within dg 
to estimate the overdensities. We find that, relative to 
the control sample, the infalling cluster members have 
~ 3 ± 2 galaxies per Mpc 2 more neighbors and contain 
~ 45 ± 7% more i?-band luminosity within dg. While 
a more accurate measurement of local density using a 
higher sampling of redshifts would be preferable, this re¬ 
sult is consistent with the idea that the infalling galaxies 
lie in higher local overdensities than the control galaxies, 
resulting in more clustering and a higher red fraction in 
the infalling sample (Figures H3l and Il4l) . _ 

The outlier rate for Q = 4 galaxies is 17% (§ IJ.2.1D . 
which may have an impact on the results presented in 
this section as a result of contamination of the cluster 
samples by field galaxies. Given that we found no red- 
shift or color dependence of the outlier rate, we expect 
both galaxies at the cluster redsliift and field galaxies to 
be affected equally by redshift inaccuracy. We use the 
CMDs to estimate the relative numbers of cluster and 
field galaxies, and find that the field galaxies comprise 
« 30% of the total at a given redshift. Therefore, we 
expect a contamination of 17% x 30% = 5% field galax¬ 
ies in the analyses above, and do not expect this to be 
sig nificant enou gh to alter our main results. 

IDressler et all (j2013i ) found that quiescent and post- 
starburst (PSB) galaxies are preferentially found in 
denser environments, including infalling groups in the 
outskirts of rich clusters at z ~ 0.4. Furthermore, they 
identified a positive correlation between the fraction of 
quiescent and PSB galaxies in the infalling groups with 






















PREPROCESSING AMONG EDISCS CLUSTER GALAXIES 


17 


increasing group mass, which they interpreted as evi¬ 
dence for “preprocessing”. While we are unable to iden¬ 
tify PSB galaxies with the LDP resolution, a prediction 
based on these results is that our infall regions contain a 
higher number of these galaxies than the field. If some 
fraction of our blue galaxies are PSB galaxies, then the 
quiescent plus PSB fraction in the infall regions would 
be even more different than the field value. Higher reso¬ 
lution spectroscopy is needed to test this hypothesis. 

In addition to environment, star fo rmat ion is corre¬ 
lated wi th galaxy stel lar mass i e.g. iKauffmann et all 
2Of). 'll iPasauali et all 120091 : iPeng et all I2010D . 
Thomas et all (|2010f) find that stellar mass may 
even be entirely responsible for the formation of early- 
type galaxies, at least at masses >10 n M 0 . While the 
galaxy stellar mass function of clusters is not correlated 
with environment (i.e^ either with in the virial region, 
outskirts, or field; iVulcani et, all I2013TI . we explore 
whether the clustered and isolated galaxies (defined 
by a separation dg) have different mass distributions. 
Following the prescription of I Bell et ah] (l2003l f. we 
estimate stellar masses using 5-band mass-to-light 
ratios, (M*/L)b, that are derived from B — V rest-frame 
colors via 

log(M*/L) B = 1.737(5 - V") - 0.942, (7) 

after converting the magnitudes to the Vega system. This 
as sumes a diet Sa l peter initial mass function as defined 
in lBell fc de .Tonal (12001h . Using Mb = 5.45 for the Sun, 
a galaxy with Mb = —19.5 and B — V = 1 has a stellar 
mass of log(M*)[M q] = 10.8. Since our infalling galax¬ 
ies span a mass range from log(M^)[M 0 1 = 9-12, we 
repeat the analyses of Figures [Td] and [fi] but restrict¬ 
ing to a narrower mass range from log(M*)[M 0 ] = 10- 
11, which is roughly symmetric about the median mass 
(log(M*)[M 0 ] ~ 10.5). For this narrower mass range, 
infalling red galaxies are significantly more likely to have 
a neighbor than blue galaxies, by « 12 ±5% compared to 
« 2±2% for the control sample. However, we cannot sta¬ 
tistically conclude that infalling galaxies of either color 
are more likely to have a neighbor of either color than the 
control sample, finding A F(< dg) « 2 ± 2%. As when 
considering the full mass range, the infalling galaxies 
with a neighbor have a higher red fraction than isolated 
galaxies, by « 50 ± 12% compared to ~ 10 ± 4% for the 
control sample. An even narrower choice of masses than 
this leaves the results qualitatively unchanged, although 
the number of galaxies becomes too few to reach sta¬ 
tistically significant conclusions like those listed above. 
While larger numbers of galaxies would help conclusively 
rule out a significant mass effect, based on these results 
we conclude that the enhanced red fraction among clus¬ 
tered galaxies is consistent with a primarily preprocessed 
origin. 

Until now we have assumed that all galaxies in the 
infall regions are falling in for the first time. How¬ 
ever, at these cluster-centric radii there exists a pop¬ 
ulation of galaxies that have already passed through 
the virial region in the past, so-called “backsplash” 
galaxies (jGill et al.i 120051 ). These galaxies may have 
been quenched on their initial passage (or passages) 
through the main body of the cluster, independent of 
any prepr oce s sing, a nd therefore must be accounted for. 
iBalogh et al.l (|2000f ) suggest that as many as 54 ± 20% 


of galaxies at distances between R 200 and 2 i ?200 are 
members of this backsplash population. Other stud¬ 
ies involving backsplash galaxies focus on distances be¬ 
tween the virial radius and 2.5 times the viri al ra¬ 
dius le.g.. iMamon et all 12004 lOman et all 120131 ) . N- 
body simulations show satellite galaxies ejected from the 
host halo cou l d com prise ~ 10% of galaxies at 2-5i?200 
(jWetzel et al.l I2013T) . While these studies show that a 
sizable fraction of backsplash galaxies could be present 
in the infall regions, these are not be expected to signifi¬ 
cantly clustered. This is at odds with our findings in Fig¬ 
ure flSlandfldl in which the red galaxies are significantly 
more clustered, and the red fraction of clustered galaxies 
is significantly higher than those without close neighbors. 
However, we cannot quantify the significance of the back¬ 
splash population in driving the enhanced quenching we 
interpret in the infall regions. It will be quite difficult to 
disentangle these two populations. 

6. CONCLUSION 

We present a spectroscopic survey of 21 EDisCS clus¬ 
ters at 0A < z < 0.8 using LDP/IMACS low-resolution 
spectroscopy. This survey contains 32,277 galaxies (with 
1, 602 within ±0.02 of the corresponding cluster redshift) 
and a precision of a z = 0.006. 

We have isolated the galaxies in the infall regions of 
these clusters using the LDP data and a simple model of 
secondary infall. The projected distance which encom¬ 
passes the infalling galaxy population, i?i n faii, agrees to 
simulations within ~ 30%. The predicted cluster ve- 
locity dispersions at z = 0 agree with the models of 
iPoggianti et al.l (|2006l) to 10%. 

With the LDP data, we identified the number of galax¬ 
ies in the infall regions, and estimate that ~ 30-70% of 
the z = 0 cluster population lies outside the virial radius 
at z ~ 0.6, a result that is not sensitive to the mass of the 
cluster over the range of cluster mass investigated here. 
This result demonstrates that studying the infalling pop¬ 
ulation is crucial to understanding how a significant por¬ 
tion of the galaxy population evolves. Furthermore, the 
ratio of the number of infalling galaxies to cluster galax¬ 
ies is typically ~ 0.3-1.5. The full range of this ratio 
is ~ 10-300%, highlighting the large cluster-to-cluster 
variation that exists. 

The red fraction in the infall regions is intermedi¬ 
ate to that in the field and clusters, suggesting that 
the process of quenching star formation has begun out¬ 
side th e virial radius, an effect previously m easured at 
z ~ 0 (|Lewis et all 120021 : iGomez et al.l 120031) . Further¬ 
more, galaxies in the infall regions show enhanced clus¬ 
tering, and the more highly clustered galaxies also show 
an elevated red fraction. These trends are indicative of 
“preprocessing”, in which galaxy star formation is shut 
off in local galaxy overdensities prior to the incorpora¬ 
tion of the galaxies into the cluster, although backsplash 
galaxies may play a role. Our sample lies at 2 ~ 0.6, 
before the epoch at which significant num bers of S Os be¬ 
gin to populate the cores of clusters (jDressler et al.lH997fc 
iFasano et al.ll2000t ). so it is plausible that the SOs in those 
cores are the remnants of quenched infalling galaxies that 
we see clustered in the infall regions. Higher resolution 
imaging is required to identify the morphologies of these 
possible progenitors. 

This dataset enables the direct study of galaxies in the 
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Fig. 13. — {Top panel) The fraction of infalling red and blue galaxies (red and blue circles, respectively) with at least one other infalling 
red/blue galaxy within a projected distance dg. Also shown are the fractions for a control sample consisting of galaxies at the same redshifts 
and radial distances as the infall regions (red and blue squares). The control data are slightly offset in dg for clarity only. (Bottom panel) 
Similar to above, only showing the fraction of infalling and control galaxies having a neighbor of either color within a distance do. 
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Fig. 14. — The red fraction of galaxies with or without a neighbor within a projected distance dg, for both the infalling and control 
samples. Galaxies with close neighbors are more likely to be red in both samples, but in the infall region the enhancement in red fraction 
is higher than for the control sample. 
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infalling regions of moderate mass clusters at interme¬ 
diate redshifts. Future work on a variety of questions 
in galaxy evolution, such as the frequ ency o f infalling 
galaxies brighter that the cluster BCG (Must et al.ll20l4 
in prep) is now possible. In particular, further studies 
seeking to understand the mechanisms that halt star for¬ 
mation in dense environments should target not just the 
virialized regions of clusters, but the outskirts as well. 
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